Dynamics and Berry Phase of Two-Species Bose-Einstein 

Condensates 

Z.-D. Chen 1 - 2 - 3 , J.-Q. Liang 1 , S.-Q. Shen 4 , W.-F. Xie 2 
1. Institute of Theoretical Physics and Department of Physics, Shanxi 
University, Taiyuan, Shanxi, 030006, China 
2. Department of Physics, Guangzhou University, Guangzhou, Guangdung, 

510405, China 

3. Department of Physics, Jinan University, Guangzhou, Guangdung, 510632, 

China 

4- Department of Physics, The University of Hong Kong, Hong Kong, China 
In terms of exact solutions of the time-dependent Schrodinger equation 

for an effective giant spin modeled from a coupled two-mode Bose-Einstein 
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between two hyperfine states of the BEC we obtain analytic time-evolution 
formulas of the population imbalance and relative phase between two compo- 
nents with various initial states especially the SU(2) coherent state. We find 
the Berry phase depending on the number parity of atoms, and particle num- 
ber dependence of the collapse-revival of population-imbalance oscillation. It 
is shown that self-trapping and phase locking can be achieved from initial 
SU(2) coherent states with proper parameters. 
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I. INTRODUCTION 

The experimental discovery of Bose-Einstein condensation in trapped atomic clouds 
opened up the exploration of quantum mechanics of mesoscopic systems in qualitatively 
new regime. The cold gas clouds have many advantages for investigations of quantum 
phenomena and hence become a test ground of quantum mechanical principles as well as 
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the interplay between macroscopics and quantum coherence. The observation of matter- 
wave interference implies the realization of coherent atomic beams, atomic Josephson effect 
and a variety of quantum interference phenomena [1]. In particular, recent experiments 
on two-component Bose-Einstein condensates (BECs) in 87 Rb atoms [1,2] have stimulated 
considerable interests in the phase dynamics and number fluctuations of the condensates. 
Aside from its intrinsic appeal, the capability demonstrated by the recent experiments might 
lead to applications also on quantum computation. Based on the macroscopic wave func- 
tion approach it is demonstrated that the Josephson effect exists in a driven two-state 
single-particle BEC in a single trap. The macroscopic quantum self-trapping as well as the 
7r-phase oscillations in which the time averaged value of the phase difference is equal to n 
have been studied extensively [3-7]. It is also shown that the population oscillation is mod- 
ulated by the collapses and revivals due to the quantum nature of the system [8-14]. The 
relative phase of two condensates in different hyperfine atomic states can be measured [2] 
using Ramsey's method of separated oscillating fields [15] and it is evident that the phase 
locking indeed occurs for small separation between condensates [1], implying the broken 
gauge symmetry. Most theoretical studies are focused on semiclassical analysis and a full 
quantum mechanical formulation of the dynamics of the two-component BECs coupled by 
time-dependent driving is certainly of interest and importance. It is well known that the 
system of two-component BECs can be described by a giant or mesoscopic pseudo-spin [16] 
using Schwinger realization of angular momentum operators in terms of two-mode bosons. 
Berry phase [17] emerges naturally in the mesoscopic pseudo-spin model [16] if the coupling 
between two components varies with time cyclically and adiabatically. The Berry phase in 
the mesoscopic spin model for the coupled two-component BECs has been explored recently 
in an elegant way by means of geometric evolutions [16]. Trapped atomic BECs make it 
possible to create mesoscopic quantum objects containing of the order 10 6 atoms in the 
same quantum state with a longer life time allowing the implementation of adiabatic evolu- 
tion which is required for the Berry Phase. However, the dynamics of the mesoscopic spin 
modeled from the two-mode BEC has not yet been studied in the quantum mechanical for- 
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malism. We in the present paper use the exact solution of the time-dependent Schrddinger 
equation for the mesoscopic spin to provide a quantum mechanical evaluation of the phase 
dynamics and the number fluctuation. With the time-evolution operator obtained by means 
of the generator of time-dependent SU(2) coherent states [18] we are able to derive analytic 
time-evolution formulas of both the population imbalance and relative phase between the 
two-component BECs for various initial states, in particular, the SU(2) coherent state with 
which the new effect of particle number dependence is discovered. Moreover, our approach 
has advantage to obtain the phase dynamics and number fluctuation for both cases with 
and without the nonlinear interatomic collisions in the same framework so that the effects 
of interatomic collision can be recognized explicitly by comparing the results between two 
cases. We show that interatomic collisions do not affect the Berry phase but lead to the 
damping and collapse-revival of the population-imbalance oscillation depending explicitly 
on the coupling strength and the total number of atoms as well. The SU(2) coherent states 
are the most realistic initial states for the two-species BEC created by coupling two hyper- 
fine states of atoms with radiation field [12]. To our knowledge we in this paper report, 
for the first time, a full quantum mechanical evaluation of the dynamics of the two-species 
BEC described by an explicitly time-dependent Hamiltonian with the initial SU(2) coherent 
states and the novel phenomena such as self-trapping, phase locking, and collapse-revival 
are recovered theoretically in the same formalism. 

The plan of the paper is as follows. In Sec. II we derive the mesoscopic spin model from 
the two-component BECs. The SU(2) coherent states, which are considered as most practical 
initial states for the two-species BEC, are briefly introduced. The dynamics and Berry phase 
of the pseudo-spin are investigated in terms of the SU(2) coherent state technique for both 
cases with and without the interatomic collisions in Sec. III. 
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II. MODEL AND INITIAL STATE 



We consider Bose-Einstein condensate of trapped atomic gas in a single trap consisting 
of two internal states which are coupled by a spatially uniform radiation field with a Rabi 
frequency A. Atoms are subjected to trapping potential V; (l=a,b). The atoms interact via 
elastic two-body collisions with the interaction potential of 5-function type. In the formalism 
of second quantization, the system is described by the Hamilton operator 

H=Y.H l + H int + H f (1) 

l=a,b 



% = J rf 3 r|§+(r) 



f ,(r) + |*+(r)*+(r)f ,(r)*,(r) | (2) 



%nt = q a , b I d 3 r^(r)§+(r)§ a (r)§ 6 (r) (3) 

H f = A(t) J d 3 r ($+(r)$ 6 (r)e^W + $+(r)$ a (r) e -^)) (4) 

We here have used the field interaction representation in rotating frame. The Rabi frequency 
A(t) is time-dependent in the sense that it can be turned on and off adiabatically [3]. The 
phase <f(t) due to the small detuning of external field from resonance excitation varies with 
time slowly and therefore we work on the adiabatically time- varying Hamiltonian. The phase 
<f(t) which we will see plays a central role in generating of the Berry phase. 

In the two-mode approximation of condensation such that $ a ( r ) ~ 2</> a (r), $&(r) ~ 60&(r) 
where a, b are the annihilation operators obeying the usual boson commutation relations, 
we have 

H a = uj a a + a + —a + a + aa (5) 
H b = uo b b + b + ^b + b + U (6) 

with 



[t), I — a,b 



(7) 



and 



Vi = qi (i 3 r|^(r)| 4 , I = a,b 



The interaction operator between two species of atoms is 



H^t = xa + ab + b , 



(8) 



where 



X = Qa,b J d 3 r\<p a (r) 



\ 2 \Mr)\ 2 . 



The transition operator induced by the external field is 



H f = G(t)(a+be ilp(t) +b + ae~ iip(t) ) 



(9) 



with 



G(t) = A(t) J d 3 r^ a (r)Ur). 



The time-dependent coupling drive is characterized by its Rabi frequency A(t) and the phase 
<p(t). Here we consider the adiabatically varying phase, <p(t) ,such that its time-derivative 
is negligibly small. The study of dynamics of the system would be greatly simplified by 
introducing of the pseudo-angular momentum operators in terms of Schwinger relation, 



J x = ^{a + b + b + a) 



(10) 



J y = ~(d + b -b+a) 



(11) 



J z = -(a + a-b+b) 



(12) 



The Casimir invariant is 



where N = a + a + b + b is the total number operator, which is a conserved quantity and thus 
is set equal to the total number of atoms N=2j with j being the quantum number of angular 
momentum. The Hamilton operator apart from a trivial constant reads 

H = u J z + qJ 2 z + G{t){J + e 1 ^ + J.e-**®) (14) 

where u = u a - uj b + (N - l) 2 ^, q=^^ - \ , and J± = J X ± U y . 

The relative phase of two-mode BEC surely can be abstracted from the expectation value 
of the angular momentum operators J± which along with the expectation value of J z (giving 
rise to the population imbalance between two components of BECs) is measurable quantity 
in experiment. 

Since the two-component BECs are experimentally created by coupling two hyperfine 
states with radiation field, it has been shown that the prepared initial state can be a par- 
ticular case of the SU(2) coherent state [12] (or known as atomic coherent state in quantum 
optics) which describes a state with a well defined relative phase between the two species. 
However, a full quantum evaluation of the dynamics for the two-species BECs coupled by 
time-dependent driving with the initial spin coherent states has not yet been given. SU(2) 
coherent state is defined as 

J • n|n) = j|n) (15) 

where n =(sin 6 cos <fi, sin 9 sin 0, cos 6*) is a unit vector. The SU(2) coherent state can be 
generated from an extreme Dicke state such that 

si(e,4>)\j,j) = \n) (16) 

where 

n{e,4>) = e^ ? -^- J+e '^ (17) 
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This is called the coherent state in the north pole gauge, as compared to the generation of 
the coherent state from the extreme state \j, —j) where it is called the state in the south 
pole gauge. The Dicke states are defined as usual J z \j,m) = m\j,m) and can be generated 
from the vacuum by boson creation operators, i.e. 

\j,m) = 1 (a + y +m (b + y- m \0) (18) 

yu + m )-U - m Y- 

The coherent state of eq.(16) can be expanded in terms of the Dicke states, 

i 

In) = £ (. ^) 2 (cos^)^ m (sin^y- m e^-^|j,m) (19) 
It is easy to verify that 



(J z ) = (n|J» = ^ cosfl, (J + ) = ysin^, (J_) = y shitfe^ (20) 

and the phase e 1 ^ is seen to be the relative phase of the two species prepared in the initial 
SU(2) coherent state. In this paper the generator of SU(2) coherent states eq.(17) is used as 
a unitary transformation to formulate the dynamics of the mesoscopic spin system modeled 
from the two-component BECs. 



III. DYNAMICS AND BERRY PHASE 
A. The case of q=0 

We first of all consider the case of q = achieved by the condition, ria + rib = % ; which, 
although a special case, is practical in the range of BEC parameters. The equality such 
that q a ~ qi, ~ q a b can be fulfilled practically for the two-component BECs consisting of 
87 Rb atoms with different internal states since the scattering lengths of atoms with the two 
internal states are known at the 1% level to be [2] in proportion a a : a a b : — 1.03 : 1 : 0.97 
where a a and cib are the same-species scattering lenghts and a a b is the scattering length for 
interspecies collisions. The ground state wave function is same, </> a (r) = <pb (r), and therefore 
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the condition rj a ks rj b x can be satisfied. The model in the case of q = is exactly solvable. 
We start with a generalized gauge transformation [19] in terms of the time-dependent unitary 
transformation [20] defined by 

R(t) = etV-e-W-J+e™) (21) 

which has the same form as the generator eq.(17) of SU(2) coherent states and is the key 
point of the present formulation. The time- dependent Schrodinger equation is convariant 
under the gauge transformation [19] such that 

ij t W{t))=H'W{t)) (22) 

where 

H' = RHR + - iRjrR + , |V>'(*)> = %W> ( 23 ) 
and the state \ip(t)) is assumed to be the solution of original Schrodinger equation 

<^|^)> = ^IV'(*)> (24) 

The auxiliary parameter A ,which is time-dependent in general, is to be determined by 
requiring that the Hamilton operator H' is diagonal in the J z representation. Using the 
relations given in appendix (eqs.Al-A4) [20], and noticing the adiabatic condition that 
^ ~ and ^ ~ 0, we obtain the Hamilton operator 



H' = a{t)J z a(t) = ^ujI + AG 2 (t) (25) 

with auxiliary parameter A chosen as 

smA = ^cosA, (26) 

UJq 

and 

cosA=^r. (27) 
a(t) y ' 



It is easy to obtain the exact general solution of the original Schrodinger equation 

\m)= E c m e- ia -U\j,m(t)), \j,m(t)} = R + (t)\j,m) (28) 

m=-j 

with \j,m) being the usual eigenstate of angular momentum J z that J z \j,m) = m\j,m). 
The total phase is seen to be 

OL m {t) = e m (t) + 7 m (£) = m \ a(t')dt (29) 

Jo 

which consists of the dynamical part given by 

e m (t) = f\j,m(t')\H\j,m(t'))dt' (30) 
Jo 

and the geometric part i.e. the Berry phase 

lm{t) = -i j\j,m(t')\-^\j,m(t')}dt' (31) 

which is defined in the usual way. In the following we only consider a time-independent Rabi 
frequency G for the sake of simplicity. For a variation of one period T i.e. tp(T) — ip(0) = 2ir 
the Berry phase is found as 

7 m (T) = — m <p (1 — cos X)d(f = m— —2n (32) 

J a 

which has an obvious, geometric meaning from the viewpoint of differential geometry that 
the one form dip is exact but not closed. Where a = \JoJq + AG 2 which is a time-independent 
parameter. The Berry phase does not depend on explicit form of the function ip(t) and is 
simply m times of a solid angle with the polar angle A in agreement with the recent result 
reported in Ref. [16] in which the Berry phase is evaluated in terms of geometric evolutions for 
the coupled two-component BECs. We, however, following the original procedure of Berry 
[17] obtain the Berry phase and the exact wave function as well by solving the time-dependent 
Schrodinger equation. In our approach the time-evolution of both population imbalance 
and relative phase between two components of BECs can be investigated analytically. The 
explicit dependence of Berry phase on the parameters of two-species BECs is also given with 
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our procedure and the novel properties of the Berry phase can be explored. The geometric 
phase is actually the same as obtained in the context of SU(2) coherent state path integrals 
[21]. To study the dynamics of the mesoscopic spin it is useful to derive the explicit time- 
evolution operator such that 

\m) = u(t)\m) (33) 

where the time-evolution operator is found from the exact general solution of the time- 
dependent Schrodinger equation \ip(t)) in eq.(28) as [20] 

U(t) = R + (t)e- iatT 'R(0) (34) 

For a given initial state |"0(O)) of the system the population imbalance between two compo- 
nents can be evaluated by 

AN(t) = N a (t) - N b (t) = 2(^(0)| l(t)m0)) (35) 

where the time-dependent angular momentum operator J z (t) in Heisenberg picture is given 
in appendix (eq.(A5)) with the help of eqs.(Al-A4). For the sake of simplicity we have set 
the initial phase to zero i.e. (p(0) = 0. To have the phase dynamics we need also the time- 
dependent angular momentum J+(t) or J-(t) in Heisenberg picture and the explicit formula 
of J+(t) is shown in eq.(A6). The time-evolution of population imbalance and expectation 
value of J+(t) given by 

are evaluated for various initial states as follows. 

(l)We consider, first of all, the initial state \ip(0))i = \j,m) and obtain in terms of 
eqs.(A5), (A6) 

AN^t) = ^K 2 + AG 2 cos(at + ip(t))] (37) 

and 
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(J+)i = msinA[-cosAe-^ (t) + cos 2 \e icd - sin 2 ^ e -^^\ ( 38 ) 

Z Z 

The population imbalance exhibits a simple oscillation. It is interesting to see a fact that 
the state of vanishing imbalance can be achieved for the case of even number of particles 
(i.e. j = y is integer) with the initial state of m — 0, but not for the case of odd number of 
particles where j is half-integer and the state of m = does not exist. 

(2)For a general SU(2) coherent state |^(0)) 2 = |n) , the population imbalance is found 

as 

AN 2 (t) = —{[cos 2 A + sin 2 A cos(at + ip(t))] cos 9 — sin A cos A sin 9 cos (39) 

Z 

+ [cos 2 — cos(at + (p(t) + <f>) — sin 2 — cos(o;t + ip{t) — (j))} sin A sin 9} 

Z Z 

The self trapping with non-vanishing population imbalance takes place for the initial state 
with = and 9 = — A. The population imbalance thus reduces to 

AN 2 (t; = 0,9 = -A) = ^ cos A = ^ (40) 

2 la 

The asymmetric trap potential i.e. non-vanishing uj q is the necessary condition to achieve 
the self-trapping. The expectation value of angular momentum operator J+(t) for the initial 
SU(2) coherent state is seen to be 

( J + ) 2 = ^{ S in A[cos 2 \e icd - sin 2 L~^+Mt)} _ cos Ae -*vW] cos Q (41) 

Z Z Z 

+ [cos 4 \e Mt + sin 4 L~^+Mt)] + * sin 2 X e~^]e^sm9 
z z z 

+ [_ C os 2 ^ sin 2 ^-{e iat + e^ at+2 ^) + \ sin 2 Ae^e"* sin 9} 
z z z 

and reduces to 

( J + > 2 (0 = 0,9 = -\) = N-e-W) (42) 

a 

for the initial state with = and = — A indicating obviously the phase locking. We see 
that self-trapping of both population imbalance and relative phase of the two-component 
BECs can be obtained simultaneously from the SU(2) coherent state. For the particular 
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case of symmetric trap potential u a = u>b ( ujq = ) we have A = — | seen from eq.(27), and 
the Berry phase then reduces to 

7 m (T, ojq = 0) = m <j> d(f = m2ir. (43) 

The population imbalance is 

AN 1 (t,uo = 0) = mcos{2Gt + (p(t))] (44) 

and 

(J+M^o = 0) = m[\^ Gt - Ie-^W +2Gt ]] (45) 

The population imbalance for the initial SU(2) coherent state vanishes seen obviously from 
eq.(40) in the case of ujq = ,while 

( J + ) 2 (0 = 0,6= -A, ^ = 0) = ^e"^) (46) 

The relative phase of two components is locked exactly to the phase of external field. The 
phase locking, which we will see, remains in the case with interatom collisions i.e. the non- 
vanishing q. It may be worthwhile to emphasize that the Berry phase of eq.(43) is trivial 
in the case of even number of particles (m is integer) while it would lead to an antiperiodic 
wave function under 2n evolution of the phase angle (p for the odd number of particles (m 
is half-integer in this case) similar to the spin parity effect in the macroscopic quantum 
coherence in spin systems [22-24]. 

B. With Nonlinear Interactions 

We now consider the general case with non-vanishing but small q for the practical model 
at hand i.e. the two-species BEC created by coupling two hyperfine states of S7 Rb atoms. For 
a single trap a reasonable condition is u = (A = — | ). With the help of time-dependent 
unitary transformation eq.(21) and the relation given by 

12 



RJlR + = T z cos 2 A + ^[(J 2 J+ + J+J z )e Mt) + (J 2 J_ + L J 2 )e^ w ] cos A sin A (47) 

+ 1 ( j2 + e i2 V (t) + j + j_ + j_ j + + J2 e -^ W) gin 2 A 

we obtain apart from a trivial constant 

S» = -GJ,-|^ (48) 

where the two-photon transition terms (proportional to J^_e t2ip ^ and J^e~ t2lp ^ ) have been 
neglected as a reasonable approximation [13,14] which is good enough for the small q in 
comparing with the transition coupling between two components, namely, q <C G. The 
Berry phase in this case is the same as eq.(43) due to u = 0. The interatom collisions do 
not affect the Berry phase. This observation for the Berry phase is new at least for the 
model we considered. The time-evolution operator is 

U q (~, t) = R + (-^ t)e« aT - + iftR(-^ 0) (49) 

The population imbalance for the initial state \tpi(0)) = \j,j) is able to be evaluated with 
the help of eq.(16) and the Dicke-state representation of the SU(2) coherent state eq.(19) as 
——l 

A Nl (q,t) = C-f £ (- 1 ) N ~ 2m L + N m + \j + m + l) l >cos{[G + q(m+ 1 -)]t-< P (t)} 

(50) 

To compare the time-evolution of population imbalance obtained here with that in the case 
of q = i.e. eq.(44) (for m — j — y ) where the time variation of imbalance is a simple 
oscillation, the time-evolution of eq. (50) is shown in Fig. 1 with various values of the ratio ^ 
and the number of particles N. It is seen that the damping goes faster when the number of 
particles, , and the ratio ^ increase (Fig.l (a),(b)). The A^-dependence of the oscillation 
of population imbalance for fixed ratio ^ is shown in Fig. 1(a) for A^ = 10 2 (dash line), 
A^ = 10 3 (dotted line), and Af = 10 4 (solid line) respectively. The nonlinear interaction 
dependence for fixed number of particles, A^ = 10 3 , is shown in Fig. 1(b) for ^ = 0.01 ~ 0.1. 



13 



With the similar method the time-evolution of population imbalance from the initial spin 
coherent state eq.(16), | ^2 (0) ) = |n), with nonvanishing q is obtained as 

AN 2 (q,t) = - ]T 

m=— ■ 

x cos{[G + qim + ^)}t - (f) - <p(t)} (51) 

where we again use the Dicke-state representation of SU(2) coherent state eq.(19). The 
oscillation of population imbalance is shown in Fig. 2 with various values of the coupling 
strength ^ and angle 9 (for the sake of simplicity we set = (p(t) = 0). Besides the 
damping the most important effect of the nonlinear interaction with initial spin coherent 
state is the collapse-revival of the population- imbalance oscillation. The particle number 
dependence of collapse-revival is shown in Fig. 2 (a) for fixed ^ and the parameter 9. We 
observe an interesting phenomena that the frequency of the collapse-revival depends on both 
the number of particles and the coupling strength ^. In Fig. 2 (b,c) we show the frequency 
behavior of the collapse-revival varying with the product of particle number and coupling 
strength , N^k. The frequency is almost the same for the same value of the product 
while different individual values of N and ^. Moreover iV^-dependence of collapse- revival 
frequency is not monotonic. A critical value N^k = 65 is found at which the frequency of 
collapse-revival approaches a minimum. To see the effect of nonlinear interatom collisions 
closely we look at the population imbalance eq.(39) with <j) = ip{t) = for the case of 
q = as a comparison. In that case the population imbalance of eq.(39) reduces to a simple 
oscillation such that 

N 

AN 2 (t,u = 0) = — cos0cos(2Gt) (52) 

for ujq = 0. It is obviously that the nonlinear interaction results in both damping and 
collapse-revival of population-imbalance oscillation. The simple oscillation of eq.(52) with 
Rabi frequency is in agreement with the experimental observation [2]. 

Particularly the initial state can be prepared such that 9 — | , the population imbalance 
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U + m + l) ( T + ni+1) * cos " +2m+1 (^) sin— 



vanishes for both the cases with and without the nonlinear interatom collisions (see eq.(52)). 
While the phase locking state is achieved again by seeing that 

(J4q,tM4> = 0,e = ^) = ^e-*M (53) 

which is the same as eq.(46) for the case of q — 0. In other words the phase locking is 
independent of the nonlinear interaction. The expectation value, (J+(q, t)) 2 , for general <fi, 
is also derived analytically with the help of eqs. (16), (19). The resulting formula is tedious 
and may not be of interest to be presented here. 

IV. CONCLUSION 

Using the exact solution of time-dependent Schrodinger equation the population imbal- 
ance and phase dynamics are evaluated with various initial states particularly with SU(2) 
coherent states which are the most realistic initial states for the two-species BEC created 
by coupling two hyperfme states of atoms with radiation field [12]. We conclude that the 
self-trapping can be achieved from initial SU(2) coherent state with asymmetric trap po- 
tential only. The phase locking is obtained also from the initial SU(2) coherent state and 
is independent of the nonlinear interaction, which may be observed experimentally in terms 
of Ramsey's method measuring the relative phase of two components of BEC in different 
hyperfme atomic states. The nonlinear interatom coupling does not affect the Berry phase 
but leads to the damping and collapse-revival of the population imbalance oscillations. The 
interesting particle number dependence of the collapse-revival and Berry phase as well is 
explored. 

Acknowledgment: This work was supported by the Natural Science Foundation of 
China under Grant Nos. 10075032. 

Appendix: 

It is easy to prove the following useful relations [20], 
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RJ + R + = J + cos 2 ^ - Le' 2 *® sin 2 ^ - J z e- 1 ^ sin A (Al) 

= J_ cos 2 ^ - J+e 2 ^ sin 2 £ - J 2 e^ w sin A (A2) 

Z Z 

RJ Z R + = J z cos A + ^(J+e 1 ^ + Le~ 1 ^) sin A (A3) 

iR^-R + = 2 d 4- sin 2 -J z --^f sin A(e*J + + e"* J_) + -^(e^ J+ - e^JL). (A4) 

St (it 2 2 dt y 1 2dt y + ' y ' 

The time-dependent angular momentum operators J z and J + are obtained respectively 

as 



J z (t) = U + (t)J z U(t) = [cos 2 A + sin 2 Acos(at + <p(t))]J g (A5) 

+1 sin A [- cos A + cos 2 -e-^W) - sin 2 -e^ + ^»] J + 
2 L 2 2 J + 

+1 sin A [- cos A + cos 2 £ e *(°*+¥>(*)) - sin 2 ^e"^^] J_, 

_ Z Z 



and 



J+(t) = sin A[cos 2 ^e iaW - sin 2 ^ e - l[aW+2 ^ )] - cos \e^ iLp{t ^]J z (A6) 

+ [ cos 4 ^ e -(t) + sin 4 ± e -i[«{t)+Mt)] + I sin 2 Ae -^)i j 
1 2 2 2 J + 

+ {_ cos 2 * sin 2 ^Mt) + e -i[a(t)+2 V ( t )]] + I ^2 Aer M*)}J_. 

2 2 2 
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Figure Captions 

Fig.l: Time-evolution of population imbalance from initial state |"0i(O)) = ip(t) = 

0. (a)§ = 0.01; N = 10 2 , 10 3 , 10 4 . (b)N = 10 3 , § = 0.01 ~ 0.1. 

Fig. 2: Time-evolution of population imbalance from initial state ^(O)) = |n), <p(t) = 
= 0, # = 0.8f . (a)§ = 0.01, iV = 10 3 ~ 8 x 10 3 . (b) A^^ dependence of AN 2 (q,t), for 
N± = 60 ~ 70, Left N = 2 x 10 3 ; Right N = 3 x 10 3 . (c)A^^ dependence of AN 2 (q, t), for 
ATi = 60 ~ 70, Left N = 4 x 10 3 ; Right JV = 5 x 10 3 . 
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